Temperature and coupling dependence of the universal contact intensity 

for an ultracold Fermi gas 
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Physical properties of an ultracold Fermi gas in the temperature-coupling phase diagram can 
be characterized by the contact intensity C, which enters the pair-correlation function at short 
distances and describes how the two-body problem merges into its surrounding. We show that the 
local order established by pairing fluctuations about the critical temperature Tc of the superfluid 
transition considerably enhances the contact C in a temperature range where pseudogap phenomena 
are maximal. Our ab initio results for C in a trap compare well with recently available experimental 
data over a wide coupling range. An analysis is also provided for the effects of trap averaging on C. 

PACS numbers: 03.75.Ss,03.75.Hh,74.40.-n,74.20.-z 



The "contact" C, introduced by Tan 1, 2] to charac- 
terize the merging of two-body into many-body physics 
in systems like ultracold Fermi gases with a short-range 
inter-particle interaction, has attracted much interest 
lately |3i-t7i|. This is especially relevant in the context 
of the BCS-BEC crossover, whereby a smooth evolution 
occurs jointly for the two- and many-body physics, from 
the presence of Cooper pairs with an underlying Fermi 
surface in the BCS limit, to the formation of molecular 
bosons with a residual interaction in the BEC limit. 

Recently, the contact C was measured in an ultracold 
gas of trapped fermionic ('^''K) atoms f8'|, from the large- 
momentum tail of the momentum distribution as well 
as from the high-frequency tail of the radio-frequency 
signal following an earlier suggestion f9|. These mea- 
surements were done from about unitarity (where the 
scattering length ap of the Fano-Feshbach resonance di- 
verges) to deep inside the BCS region (more precisely, 
in the coupling range —2.5 ^ {kpciF)^^ ^ -1-0.5 where 
kp = {2mEFy^^ is the Fermi wave vector expressed in 
terms of the Fermi energy Ep = u;o(3A^)^/'^, m being the 
atom mass, N the total number of atoms, and wq the 
average trap frequency [lO|), and in a temperature range 
about the critical temperature T^. 

Several questions concerning the contact C remain 
open. They hinge on the recent experimental results of 
Ref. 8] (like the temperature and coupling dependence of 
C, and the effect that trap averaging has on the values of 
C), as well as on more theoretical issues. These include 
the identification of the (approximate) spatial boundary 
between short- and medium-range physics that can be 
associated with the contact C, the effects that improved 
theoretical approaches have on the values of C, and the 
interconnection with the presence of a pseudogap in the 
single-particle excitation spectrum. 

In this Letter, we address these questions and calcu- 
late the contact C using a t-matrix approach [ll'| that 
proved successful in comparison with ARPES-like data 
for ultracold Fermi atoms Il2ll. and also using a nontriv- 
ial extension of this theory jl3| which takes into account 



the residual interaction among composite bosons. This 
is to verify to what an extent improvements on the de- 
scription of the medium-range physics (over and above 
the results of the t-matrix) influence the values of C in 
different coupling and temperature ranges |14{ . 

We shall, specifically, be concerned with the temper- 
ature dependence of C over an extended temperature 
range up to (several times) the Fermi temperature Tp, 
to determine how the value of C is affected by the 
pseudogap physics extending above Tc in the unitary 
(—1 ^ {kpap)~^ < -1-1) regime, and to address the re- 
lated question of how trap averaging influences the value 
of C with respect to that of a homogeneous system with 
the same nominal temperature and coupling. 

The contact C was originally introduced to account 
for the large wave- vector behavior of the fermionic dis- 
tribution n{k) (for spin component) of the homogeneous 
system, such that n{k) « Chfc""'. Here, the suffix 
"h" stands for homogeneous, k — |k| is in units of 
kp = (STT^n)^/"^ where n is the total particle density, and 
n{k) is normalized such that Jj^^ n{k) = 1/2. Alterna- 
tively, Ch can be extracted from the high-frequency tail 
of the radio- frequency (RF) spectrum Inpiuj) per unit 
volume, so that Irp{uj) w {Ch/2^^^Tr^)oj~^/^. Here, the 
frequency uj is in units of Ep and the RF spectrum is 
normalized such that /_ duj Ijip{uj) = 1/2. The above 
asymptotic form of Ib,p(u}) holds provided final-state ef- 
fects can be neglected |9l.ll5lll6|. 

Similar asymptotic behaviors can be obtained for the 
trapped system. Preserving the above normalization as 
for the homogeneous system, we write within a local- 
density approximation; 
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Here, the suffix "t" and kp refer to the trapped system, 
the spatial position r is in units of the Thomas-Fermi 



radius Rp — [2_Ei?/(ma;Q)]^/^, and n(r) and Ch(r) are 
the density and contact locally in the trap. By a similar 
token, the large-oj behavior of the total RF spectrum of 
the trapped system reduces to: 
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with the overall normalization of the homogeneous case. 

This local-density analysis shows that, by adding con- 
tributions of different shells with a weight proportional 
to n(r)''/^, the trap modifies the value Ch for a homo- 
geneous system at the same nominal coupling {kpaF)~^ 
and (relative) temperature T/Tp. In this respect, the 
temperature dependence of Ch provides an important in- 
formation, because different shells are at different relative 
temperature with respect to the local value of Tp . 

We begin our analysis with the homogeneous case at 
Tc and show in Fig. [Ha) the coupling dependence of Ch, 
when n{k) is obtained within the t-matrix approximation 
of Ref.|lll| (full line) and its improved Popov version of 
Ref. 13] (dashed line). The calculations are done at the 
value of Tc of the respective theories. The plot also shows 
(light-dotted lines) the leading (low-temperature) ap- 
proximations for Ch obtained in the BCS weak-coupling 
and BEC strong-coupling limits, which are given by 
the expressions (4/3)fc^a|. and ATT/{kpap), respectively, 
with the crossover region —1 ^ {kpaF)~^ ^ +1 marking 
the change between these two limiting behaviors. 

Note how the difference between the full and dashed 
lines in Fig.[IJa), which originates from the activation of 
the residual interaction among the composite bosons, is 
appreciable only close to unitarity where particles cor- 
relate with each other within the inter-particle spacing 
kp . The smallness of this difference resulting from 
our calculations confirms the validity of the t-matrix ap- 
proximation for the contact Ch and for the high-energy 
scale to which Ch is associated (see below). Yet, this 
difference is relevant for the physical interpretation of 
the contact Ch as characterizing the effects of medium- 
range (many-body) physics over and above the short- 
range (two-body) physics. This interpretation is also con- 
sistent with rewriting Ch = (37r^/4)(Aoo/-Ei^)^ in terms 
of the high-energy scale Aoa introduced in Ref.[9|, such 
that Aoo = 2TT\ap\n/m embodies in weak coupling the ef- 
fects of surrounding particles through a mean-field shift 
[17j . while A^ = 47rn/(m^ai?) reflects a standard re- 
lation in strong coupling [1^ between the density and 
the gap parameter within BCS theory. For compari- 
son. Fig. (Ha) also reports the coupling dependence of 
Ch at zero temperature (dashed-dotted line) within the 
t-matrix approximation, to which the above approximate 
expressions in the BCS and BEC limits converge. 

Figure [T]Jb) shows the temperature dependence of Ch 
at unitarity over an extended temperature range within 
the t-matrix approximation. The rather slow decay of 
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FIG. 1: The contact Ch for the homogeneous case: (a) At 
Tc vs the coupling {kpap)^^ , obtained within the t-matrix 
approximation (full line) and its improved Popov version 
(dashed line). The leading approximations in weak and strong 
coupling (light-dotted lines) as well as the T — result within 
the t-matrix approximation (dashed-dotted line) are also re- 
ported, (b) At unitarity vs the temperature T (in units of the 
Fermi temperature Tp), obtained within the t-matrix (full 
line). The high-temperature approximation to this curve is 
also reported as a reference (dashed line) and extrapolated to 
T = 0. [See the text for the meaning of the inset.] 



Ch at high temperature is consistent with the expec- 
tation that Ch is not related to long-range order. The 
temperature behavior of Ch steepens up at low tempera- 
ture when entering the pseudogap region for T/Tp < 0.5 
[12| , thus evidencing the emergence of a local (medium- 
range) order which is also responsible for the pseudo- 
gap. Such an enhancement of the value of Ch appears 
most evident when the calculation is continued below Tc 
|l9| . with the result that a cusp appears in Ch at Tc 
where the effect of the pseudogap is maximum. To em- 
phasize this enhancement, we have indicated in Fig.[IJb) 
the extrapolation of the high-temperature behavior of Ch 
(dashed line) down to T = 0, on top of which the con- 
tribution associated with the pseudogap region about Tc 
appears evident. Our value (=3.23) for Ch at T = com- 
pares well with that (=3.40) extracted from the Monte- 
Carlo calculations of Ref.[2Cf|. For completeness, the 
inset of Fig. [Ijb) reports the temperature dependence 
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FIG. 2: The high-energy scale Aoo (full line) and the low- 
energy scale Apg (dashed line) are reported (in units of Ef) 
vs the coupling {kpap)'^ at Tc- 



of Ch for (kpap)^^ = —1.0 (full line/left scale) and 
(kpap)^^ = -1-1.0 (dashed line/right scale). The maxi- 
mum at T /Tp ~ 0.5 in the weak-coupling curve is consis- 
tent with the Fermi-liquid behavior discussed in Ref. 14 1 . 
The temperature dependence of Ch is here reported 
for the first time and deserves further comments. On 
physical grounds, the enhancement of Ch when enter- 
ing the fluctuative region from above Tc is due to the 
strengthening of local pairing correlations, still in the 
absence of long-range order. In the present approach, 
pairing correlations are embodied by the pair-fluctuation 
propagator in the particle-particle channel, whose wave- 
vector and frequency structures give rise to a characteris- 
tic low-energy scale Apg in the single-particle excitations 
ll|, which is referred to as the pseudogap. The con- 
tact Ch, through its alternative definition in terms of the 
high-energy scale Aoo that was previously mentioned, is 
also related to a wave-vector and frequency averaging of 
the very same structures in the pair-fluctuation propa- 
gator [9|. The inter-dependence between the two energy 
scales Aoo and Apg can then be explicitly appreciated 



in Fig. [21 where they are shown at Tc vs the coupling 
{kpap)~^. The two energy scales do not fully relate to 
each other in weak coupling where the contact is dom- 
inated by the mean-field interaction, but become very 
close in value at unitarity where strong local-pairing cor- 
relations dominate both thermodynamic (Aoo) and dy- 
namic (Apg) quantities. 

Although the values of Ch are extracted from the 
asymptotic (scale-free) power law n{k) « Ch k~^, a char- 
acteristic value fcc can nevertheless be identified at which 
n{k) has reached Ch k~'^ within, say, a few percent accu- 
racy. The length scale corresponding to k^ is approxi- 
mately the spatial range at which two-body physics inter- 
faces with medium-range physics. This range is expected 
to be about the size of the composite bosons in strong 



FIG. 3: The characteristic values of kc (in units of kp), at 
which the asymptotic power-law behavior n{k) ~ Ch k~ is 
reached within 2.5% (dashed line), 5% (dashed-dotted line), 
and 10% (full line) accuracy, are shown at Tc vs {kpap)^^ ■ 
The inset shows, correspondingly, the characteristic value ujc 
(in units of Ep), at which the asymptotic power-law behavior 
Irf{oj) ^ (Ch/2^/^7r^)a;"^/^ of the RF spectrum is reached 
within 10% accuracy. 



coupling and the inter-particle spacing in weak coupling. 
Figure |3] shows the values of kc extracted in this way 
at Tc vs {kpap)~^ within the t-matrix approximation, 
for the three distinct values (2.5,5,10)% of the above 
percent accuracy. In extreme weak coupling, our numer- 
ical values recover the limiting ones obtained for a dilute 
Fermi gas [21|. Note how fee reaches a pronounced min- 
imum for the coupling at which the chemical potential 
vanishes, away from which kc increases most markedly 
on the strong-coupling side as expected from the increas- 
ing spatial localization of composite bosons. 

In addition, the inset of Fig. [3] shows the value of wc at 
Tc vs {kpaF)~^ , extracted within a 10% accuracy from 
the large-w behavior of the RF spectrum. (Recall that 
the w^^/^ tail of the RF spectrum originates from the 
short-range behavior of the 2-body wave function |22| . 1231 1 
in the absence of final-state effects 



24|.) Comparison 



with kc with the same accuracy yields the relation 2ujc = 
k^ that holds approximately for all couplings. 

The values of the contact Ct obtained from Eq.(IT]), 
by adding the asymptotic contributions from all shells 
in the trap within the t-matrix approximation, are re- 
ported in Fig. H at Tc vs {kpap)''^ (full line). They 
are compared with the values (filled and empty circles, 
stars) obtained experimentally in Ref.fSj through alter- 
native procedures. The theoretical value obtained at 
unitarity for T = is also reported for comparison 
(empty square). Due to difficulties in extracting the 
asymptotic behavior of n{k) from experimental data, the 
figure also shows the theoretical values of Ct (dashed- 
dotted line) obtained upon averaging k'^n{k) over the 
interval fcmin l£ k < A:niax- This follows the procedure 
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particle excitations about Tc- 

Discussions with E. Cornell, D. Jin, and J. T. Stew- 
art are gratefully acknowledged. This work was partially 
supported by the Italian MIUR under Contract PRIN- 
2007 "Ultracold Atoms and Novel Quantum Phases" . 



FIG. 4: The contact Ct obtained within the t-matrix ap- 
proximation for the trapped case (full line) is shown at Tc vs 
the couphng {kpap)~^ , and compared with the experimental 
values of Ref.Q] (filled and empty circles, stars). The inset 
shows, correspondingly, Ct at unitarity vs T/Tp. 



used to extract the experimental values of Ct, for which 
^min = 1-55 when (/cfOf)"^ < —0.5 and fcmin — 1-85 
when —0.5 < {kpaF)~^, while fcmax = 2.5 (^l- 

It is worth noting that the data of Ref.[8| which are 
reported in Fig. 21 have been originally compared with 
a theoretical curve from Ref. 4], where the behavior of 
Ct at r = was obtained across the unitary regime by 
interpolating known results in the BCS and BEC limits. 
That interpolation has to be contrasted with the com- 
pletely ab initio theoretical calculation at Tc reported by 
the full line in Fig. |4l which spans a wide coupling range 
just across the unitary regime. 

Finally, the inset of Fig. |4] displays the temperature 
dependence of Ct for the trapped case at unitarity within 
the t-matrix approximation (full line), with the vertical 
arrow indicating the corresponding value of Tc- Contrary 
to the result for the homogeneous case of Fig. [TJb), no 
cusp now appears in Ct at Tc for the trapped case. An 
expanded discussion about the effects that trap averaging 
has on the values of C is provided in [251] . 

In conclusion, we have presented a detailed analy- 
sis of the contact intensity C for a Fermi gas over the 
temperature-coupling phase diagram. The values of C 
have been obtained from the large-fc behavior of the mo- 
mentum distribution n{k) as well as from the large-w tail 
of the RF spectrum Iiip{uj), both for the homogeneous 
and trapped case. For the latter case, good agreement is 
obtained with recent experimental data. The effects of 
pairing fluctuations on C have been determined by the 
t-matrix approximation and its improved Popov version 
that takes into account the interaction among composite 
bosons. We have found that the values of C are strongly 
affected by the emergence of a pseudogap in the single- 
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SUPPLEMENTAL MATERIAL: 

"TEMPERATURE AND COUPLING 

DEPENDENCE OF THE UNIVERSAL 

CONTACT INTENSITY FOR AN 

ULTRACOLD FERMI GAS" 

We discuss the connection between the values of the 
contact intensity for the whole trap and for a homoge- 
neous system. 

It is interesting to compare the values of Ct obtained 
for the whole trap with its approximation obtained from 
Eq.(2) of the main text, whereby one identifies the 



shell at 



corresponding to the maximum of the ra- 



dial weight function (32/7r) r^ [Stt^ n{r)]'^^'^ /kp, and then 
takes Ch(?'max) therein outside the integral. The result of 
this procedure at Tc is reported vs the coupling lkpap)~^ 
in Fig. 5(a), where the inset shows an example of the 
shape of the radial weight function (full line/right scale) 
and of Ch(r) (dashed line/left scale) at unitarity. The 
good agreement, which results between the calculation 
for the trap (full line) and the approximation that se- 
lects the contribution of the most important shell (dashed 
line), shows to what an extent the results of Ct for the 
whole trap can be used, together with knowledge of the 
density profiles, to extract the values of Ch for the ho- 
mogeneous case. 

The same procedure can be applied to interpret the 
temperature dependence of Ct, reported in the inset of 
Fig. 4 of the main text and reproduced here in Fig. 5(b) 
for convenience (full line). In particular, it is interesting 
to understand how the trap averaging washes out the 
peak about Tc obtained for the homogeneous case (see 
Fig. 1(b) of the main text). 

To appreciate this effect, we have reported in the inset 
of Fig. 5(b) the temperature dependence at unitarity of 
the integral of the radial weight function (full line/right 
scale) and of Ch(?'max) (dashed line/left scale) from above 
to below Tc- While Ch(?'max) retains the characteris- 
tic cusp feature of the homogeneous case (with a maxi- 
mum at T = 0.8Tc), the steady increase of the integrated 
weight for decreasing temperature more than compen- 
sates for the decrease of Ch(?'max) when T < 0.8Tc, thus 
masking eventually the cusp feature in the integrated 
quantity. 

The same approximate procedure, that resulted in the 
dashed line of Fig. 5(a), can be applied to reproduce the 
temperature dependence of Ct at unitarity above Tc, be- 
cause in this case the two functions in the integral of 
Eq.(2) of the main text have a smooth behavior similar 
to that shown in the inset of Fig. 5(a). At given T below 
Tc, however, the cusp present in C\^{r) requires us to split 



it as the sum of a smooth background C\^ '{r) and of a 
peaked contribution Cj^ (r), yielding approximately: 
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(4) 



These approximate results, from above to below Tc, are 
shown by the dashed line in Fig. 5(b) (to which the sec- 
ond term on the right-hand side of Eq.(j4|) below Tc gives 
at most a 15% contribution). 
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FIG. 5: The contact Ct obtained within the t-matrix approx- 
imation for the trapped case, is shown (full line) (a) at Tc 
vs the coupling {kpap)'^ , and (b) at unitarity vs T/Tf, and 
compared with the approximation (dashed line) that relies on 
the maximum of the radial weight function. [See the text for 
the meaning of the insets.] 



